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ABSTRACT 

The 2 Micron All-Sky Survey (2MASS) TuUy-Fisher Survey (2MTF) aims to measure Tully-Fisher 
(TF) distances to all bright incHned spirals in the 2MASS Redshift Survey (2MRS). Essential to this 
project is a universal calibration of the TF relation in the 2MASS J (1.2 /um), H (1.6 /Ltm) and K (2.2 
/im) bands. We present the first bias corrected or universal TF template in these bands. We find 
that the slope of the TF relation becomes steeper as the wavelength increases being close to L oc 
in K-band and L oc v^'^ in J and H-bands. We also investigate the dependence on galaxy morphology 
showing that in all three bands the relation is steeper for later type spirals which also have a dimmer 
TF zeropoint than earlier type spirals. We correct the final relation to that for Sc galaxies. Finally 
we study the scatter from the TF relation fitting for a width dependent intrinsic scatter which is not 
found to vary significantly with wavelength. 

Subject headings: galaxies: distances and redshifts — galaxies: clusters — galaxies: fundamental 
parameters — distance scale: infrared: galaxies 
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1. INTRODUCTION 

The 2 Micron All-Sky Survey (2MASS) Tully-Fisher 
Survey (2MTF) is making use of existing high quality 
rotations widths, new HI widths and 2MASS photome- 
try to measure Tully-Fisher (TF) distances for all bright 
inchned spirals in the 2MASS Redshift Survey (2MRS; 
Huchra et al. 2005). Previous peculiar velocity surveys 
have struggled to meet their potential because of large 
errors on individual measurements, poor statistics and 
uneven sky coverage. A survey based on 2MRS will pro- 
vide both better statistics and greatly reduce the im- 
pact of the ZoA providing a qualitatively better sample 
for velocity-density field reconstructions. A significant 
amount of new HI observing (both at the Green Bank 
Telescope, and the Parkes Radio Telescope) has already 
been completed for the 2MTF project (Masters et al. in 
prep, see Masters 2007 for more details). 

The first step in constructing the 2MTF catalog is to 
have a universal calibration of the TF relation in the 
2MASS J, H and K-bands. As discussed in Masters et 
al. (2006; hereafter M06) the estimation of proper bias 
corrections for deriving TF templates is essential, not 
only if the template is to be used as a basis for unbi- 
ased distance estimates to galaxies with a wide range of 
rotational widths, but also if the template is to be used 
as a constraint in theories of galaxy formation. A bias 
corrected, or universal TF relation in the 2MASS J, H 
and K-bands does not currently exist. In this paper we 
provide such a TF template relation for use with 2MTF 
and in other applications. We also investigate the de- 
pendence of the NIR TF relation with wavelength and 
galaxy morphology and study the extent of the TF rela- 
tion. Finally we measure the intrinsic scatter of galaxies 
from the TF template relations as a function of galaxy 
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rotation width. Wc finish by providing a recipe for de- 
riving TF distances from data consisting of galaxy mag- 
nitudes, inclinations and rotational widths. We parame- 
terize quantities dependent on the Hubble constant with 
h ~ _ffo/100 km s^^ Mpc~^ where possible, and use 
Ho = 74 km s^^ Mpc^^ when we need to compare with 
Ho independent measurements {e.g. calibrations of the 
TF zeropoint from primary distance measures). 

2. THE 2MASS TF TEMPLATE SAMPLE 

The 2MTF template sample discussed here originated 
in a cross-match of the 2MASS Extended Source Cat- 
alog (XSC) with all galaxies with rot ation widths ei- 
ther from the Cornell HI digital archive (jSpringob et al. I 
l2005f ). or from the database of optical rotation curves 
(ORCs) maintained by Martha Haynes and R iccardo 
Giovanelli at Cornell (see ICatinella et al. 1 12005' for the 
most recent description of the ORCs) . From this cross- 
match we pick galaxies in the vicinities of the 31 nearby 
clusters described in M 06. In M06 the SFI-I--)- sample 
(|Springob et al. |[2007[ ) is used as a base for their related 
I-band TF template sample. This "SFI-I— I- template 
sample" was comprised of galaxies which in addition 
to rotation widths from th e above sources have I-band 
photometry as described in lHavnes et al. iri999ai Using 
the same set of clusters as in M06 will ease comparisons 
between the I-band template from M06 and the J, H 
and K-band templates to be derived here. Using the ter- 
minology of M06 (and Giovanelli et al. 1997b; hereafter 
G97b) who describe an in sample in which all galax- 
ies are clearly bona fide cluster members and a larger 
in+ sample which may include galaxies merely in the 
vicinity of the cluster we consider mostly the in+ sample 
for these clusters. 

Roughly 70% of the SFI-H-h template sample (M06) are 
found in this "2MTF template sample" . Some SFI-h-h 
galaxies are too low surface brightness in the NIR to 
be detected by 2MASS. This is especially a problem for 
the SFI-I--I- since it was designed to be dominated by 
late type spirals which tend to have a lower NIR surface 
brightness. Of the final 2MTF template sample, 65% are 
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also in the SFI++ template sample (the remaining 35% 
have available rotation widths, are present in the 2MASS 
XSC and located near one of the "template clusters" , but 
do not have I-band photometry). 

2.1. Rotation Widths 

As described above, rotation widths are taken from 
the same source as the SFI++ template sample 
(M06), namely the HI digital archive published in 
ISpringob et al. I (|2005l) supplemented by optical rotation 



curves (see e.g. iCatinella et~aL\\200^ . The correc- 



tion from ORG to HI width (measured at 50% of the 
(peak— rms) on a fit to the sides of the profile) is de- 
scribed in M06 and is applied, with the modification 
that the optical radius (ass, which is the radius confining 
83% of the optical light) is estimated for galaxies which 
have only 2MASS photometry from the 2MASS parame- 
ter Text - the 2MASS XSC extrapolation or total radius. 
For galaxies with both 2MASS and I-band photometry 
we find as3 = 0.56rext, with a scatter of 13" - this is the 
relation used to calculate 033 for the 35% of the 2MTF 
sample which do not have available I-band photometry. 
The observed width must also be corrected for instru- 
mental effects (As), inclination, cosmological broadening 
and the effect of turbulent motions (At) in the disk: 



1 + z 



A, 
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We use At = 6.5km s ^. The instrumental correction 
will depend on the d etails of the telescope used: see 
ISpringob et~ar\ ()2005D for more information on this. 

2.2. Total Magnitudes and Inclinations 

The 2MASS XS C photometri c quant ities are described 
in great detail in ICutri et al. I (|2006l ). The quantities 
of interest here are the total magnitudes in the J, H 
and K bands (the extrapolated magnitudes are used), 
and the axial ratio. Inclinations are estimated from the 
deep I-band imaging taken as part of the SFI++ or 
SFI projects and described in .Hav nes et al. I (|1999aD for 
galaxies where this is available. In the 35% of galaxies 
for which no I-band imaging is available the J-band axial 
ratio (after correction for seeing as described in Masters, 
Giovanelli & Haynes 2003) is used to estimate the inclina- 
tion. We choose the J-band over longer 2MASS bands (or 
the co-added image) in the hope that it best traces the in- 
clination of the disk and is least affected by the bulge. We 
also use an intrisic axial ratio of 0.2 for all types of spirals 
to mitigate the effect of the bulge on the most edge-on 
galaxies. For all SFI-I-I- galaxies in the 2MASS XSC (ie. 
not just galaxies used in this paper), an approximately 
one-to-one correlation is found between the I-band and 
J-band estimates of the inclination, with a rms scatter of 
~ 8.5°. There is a slight tendency for the J-band axial 
ratio to imply the galaxy is less inclined than the I-band 
measurement, which is especially prominent for the most 
face-on (ie. lowest inclination) galaxies. We do not use 
any galaxies with J-band inclinations of less than 25° - 
galaxies with true inclinations this low would not make it 
into the SFI-I— I- as they would not have been considered 
good candidates for TF distance measurement. There 
is no significant trend of the difference between the two 
measures of inclination with galaxy morphological type. 



By looking at the scatter in the axial ratio measured in 
J-band vs. that measured in I-band we estimate that 
the error in the J-band axial ratio is ±0.1, remain- 
ing approximately flat with decreasing ellipticity. This 
is the error we carry forward into the error on the rota- 
tion widths for galaxies in which J-band inclinations are 
used. By comparison the I-band axial ratio has a mean 
reported error of ±0.04, increasing to ±0.06 for the most 
face-on galaxies, and decreasing to less than ±0.02 for 
the most inclined. The use of J-band axial ratios for 
35% of the sample will obviously add some scatter into 
the final TF template. 

Galactic extinction is estimated using the DIRBE dust 
map (Schlegel, Finkbeiner & Davis 1998), while dust ex- 
tinction internal to the galaxies themselves is estimated 
us ing the lumino s ity de pendent empirical relations found 
in lMasters et al\ (I2003D. A small K-corre ction is apphed, 
also as described in [Masters' "eTZI (p003l) . 

3. CONSTRUCTING A GLOBAL TEMPLATE 

Many schools of thought exist on the best way to con- 
struct a TF template. In one method no atte mpt is made 
to co rrect for observational biases {e.g. Pi zagno et al.l 
l2007t ). The argument is that such corrections depend on 
various assumptions and so add an unmeasurable source 
of error into the template. Such templates must not be 
considered global in any sense and can only be compared 
to theory when mock catalogs with the same selection 
criteria are constructed. They also have the distinct dis- 
advantage of introducing obvious biase s into distances 
deriv ed from TF. Another method {e.g. iTheureau et al.l 
|2007( ) attempts to construct an "unbiased" sample of 
galaxies so that no corrections are necessary - this has 
the disadvantage of requiring a large part of the data be 
rejected from the template sample, and because the re- 
jected galaxies lie preferentially at the low mass end of 
the relation this will also result in a TF template relation 
which is not completely universal since it will be fit only 
to the most massive spirals. It is possible to construct 
a bias corrected or universal TF template in which ob- 
servational sample biases are corrected using reasonable 
assumptions about the properties of the parent galaxy 
population. With the right assumptions this method is 
the only one which produces a universal or global TF 
relation which can be compared directly to theory and 
used without concern to derive TF distances. This is the 
approach we will take here. 

A detailed description of all the relevant bias correc- 
tions which need to be applied to derive a universal TF 
relation has recently been published in M06; we refer the 
reader to that paper for more details. Both in M06 and 
this paper, the method used is based on that described 
in G97b. What follows is a list of all bias corrections, 
including details of how they have been modified to ap- 
ply to 2MASS photometry and the template sample dis- 
cussed here. 

These corrections are applied after the magnitudes 
have been corrected for Galactic and internal extinction 
and the inclinations have been corrected for seeing; and 
after the HI widths have been corrected for instrumental 
broadening, the ORG widths corrected to the HI scale, 
and all widths have been corrected for inclination and 
cosmological broadening. 

1. Incompleteness bias correction: In any given clus- 
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ter we preferentially observe galaxies which are 
scattered above the TF slope (because they are 
brighter and therefore more likely to make it into 
the sample) . This has the result of shallowing the 
TF slope, brightening the zero point and reducing 
the measured scatter. The correction we apply for 
this effect is derived following the method described 
in G97b and M06. It is derived from Monte Carlo 
simulations of the difference between the expected 
and observed magnitudes requiring an input TF 
relation, TF scatter and spiral LF (which is used 
to estimate the completeness of the sample). We 
do the calculation only once, using as an input the 
TF relations and scatter uncorrected for this ef- 
fect (but already corrected for the morphological 
dependence as described in #5 below). We discuss 
below the impact of these assumptions on the size 
of the bias correction. 

A luminosity function was calcula ted for late-type 
galaxi es (T > 0) in K-band by iKochanek et al.l 
()200lD . This LF includes SO/a galaxies (T = 0) 
as well as later types, but they found that the LF 
parameters were insensitive to moving the bound- 
ary to include only Sa and later galaxies. They 
fit a Schechter function with Mk* — 5 log h = 
-22.98 ± 0.06 mags and a = -0.92 ± 0.10. This 
is the luminosity function we use to calculate the 
incompleteness of the template sample. To each 
measure incompleteness for a cluster we fit a func- 
tion of the form 

c{y) = -. \-, (2) 

The parameters of the fit to the K-band incom- 
pleteness function for each individual cluster sub- 
sample are given in Table [TJ To convert the com- 
pleteness to J and H bands we use a constant color 
estimates oi J - K = 1.04 and H - K = 0.28 
which are the average for all galaxies in the tem- 
plate sample. There are no large trends oi H — K 
color with absolute magnitude and only a slight 
trend for J — K color to be smaller for dimmer 
galaxies (which could result in a very slight over- 
estimation of the incompleteness in J-band for the 
least luminous galaxies). This J — K color is a little 
redder than that measured for a complete sample of 
2MASS galaxies (Huchra et al. in prep.) although 
the H — K color is consistent. Galaxies in the tem- 
plate sample are larger and more inclined that the 
average galaxy in 2MASS so even though this color 
is corrected for inter nal extinction as described in 
[Masters et al\ ()2003r ) we expect them to be a little 
redder than average because this mean correction 
may still be a underestimate for the brightest most 
edge-on galaxies. 

We use a bivariate fit to the full sample as our 
estimate of the universal TF relation. We do one 
round of the iteration - we use to calculate the bias 
corrections the initial relations (not corrected for 
incompleteness bias, but after the morphological 
correction discussed below is applied) which are 

Mk - 5 log /i = -22.228 - 8.919(logW^ - 2.5), 
Mh - Slog /i = -21.954 - 8.289(logM/ - 2.5), 



Mj - 51og/i = -21.198- 8.030(logl^ - 2.5). 

This can be compared with the final output in Ta- 
bic [2l The slopes here are shallower than the final 
relation by ^ 1 mag/dex in K- and J-bands and 0.7 
mag/dex in H-band. The zeropoints are brighter 
by ~ 0.2 mag in K- and J-bands and 0.1 mag in 
H-band. As shown in G97b and [Masters I (|2005l) 
changing the zeropoint and slope of the assumed 
UTF by these amounts will not change the magni- 
tude of the derived bias correction significantly. 

The assumed scatter we use is also the scatter from 
the initial relation (uncorrected for incompleteness 
bias but after correction for morphology) of 

eint,K = 0.41-1.43(logW^-2.5), 
ei„t,H = 0.43- 1.26(logVF-2.5), 
eint,j = 0.45- 1.08(logiy-2.5), 

(see Section 4 for an explanation of the derivation 
of this and the final output values). It is this as- 
sumption which has the biggest effect on the mag- 
nitude of the calculated bias correction. However, 
the scatter from the final TF relation is only 0.05 
mag smaller than this value (in all 3 bands), and 
this is not enough of a difference to significantly 
change the derived corrected for the incomplete- 
ness bias (as shown in G97b and Masters 2005). 

The magnitude of the bias correction for a given 
galaxy depends on a combination of the abso- 
lute magnitude of the galaxy and the distance to 
the cluster in which the galaxy is located (in the 
nearer clusters the sample is obviously complete 
to a fainter absolute magnitude limit). In G97b 
and M06 the bias correction is shown separately 
for galaxies in each of the clusters used in the tem- 
plate. Here we instead show bias corrections for all 
galaxies in all clusters together (Figured)). 

2. Mean distances to cluster. This small fixed bias 
for all galaxies in a given cluster arises from the 
fact that the log of the average of the distance is 
not equal to the average of the log of the distance 
(which is what we implicitly fit for when using mag- 
nitudes). This bias is easy to calculate once a line 
of sight depth for each cluster has been estimated. 
We use as a size estimate 1 Abell radius for the 
in sample, 2 Abell radii for in+. This bias is al- 
most always much less than 0.1 mag, and is given 
for each cluster in Table [T] (/3sizc)- 

3. Cluster size - incompleteness. In a given cluster we 
preferentially observe dimmer galaxies only in the 
closest parts of the cluster - creating a small bias in 
the observed magnitudes which is easily corrected 
for using the completeness we have already calcu- 
lated for each cluster and assumptions about the 
size of the cluster. This bias is very small (less 
than 0.1 mags) for most galaxies in the sample. 

4-. Edge of catalog bias and Malmquist bias. Following 
G97b and M06 we argue that both of these biases 
are negligible in the sample and therefore make no 
correction for them. 
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Fig. 1. — Incompleteness bias corrections in J, H and K-bands 
for all galaxies. Points are coded by the redshift of the clus- 
ter the galaxy belongs to. Blue dots: fduster < 3000 km s~^\ 
green squares: 3000 km s^^< "Uclustor < 7000 km s~^; red crosses: 
f cluster > 7000 km s~^. This illustrates how the size of the bias 
correction for a given galaxy depends both on the luminosity of the 
galaxy and the completeness of the individual cluster sub-sample 
the galaxy belongs to. It also illustrates that clusters at similar 
redshifts but in different parts of the sky have quite different levels 
of completeness. 

5. Morphological correction: We observe that the 
slope for earher type spirals is shallower than for 
later types and that they have a brighter zeropoint. 
We correct to an "Sc" (type T = 5) relation in K- 
band using, 




AMsa = 0.19 
AMsb = 0.01 



3.16(logVF-2.5) 
1.46(logW-2.5) 



T < 2 
3 < T < 4 



T < 2 
3 < T < 4. 
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In H-band we use, 

AMsa = 0.21 -2.92(logW^- 2.5) 
AMsb = 0.01 -1.26(logH/- 2.5) 

In J-band we use, 

AMsa = 0.23 - 2.86(logW- 2.5) 
AMsb = 0.02 -1.31(logW^- 2.5) 

These values are as measured using the data for 
each sub-sample (separated by morphological type) 
before it has been corrected for sample incom- 
pleteness bias (above) . The final fits to the sample 
separated by morphological type once the incom- 
pleteness bias correction has been applied are listed 
in Tabled! and shown for K-band in Figure^ They 
are very similar to what is used here. 

The variation of the TF relations with morpholog- 
ical type is found to have a similar magnitude in 
all 3 bands, and a similar dependence on the ro- 
tational width. The morphological dependence of 
the TF relation is discussed further in Section [3721 
The caveat to be added here is that if the incom- 
pleteness for earlier type spirals is greater than for 
the later types, then this variation in the TF rela- 
tion may not be physical but rather the impact of 
an underestimate of the incompleteness correction 
for earlier type spirals. Naively, in the NIR we ex- 
pect the selection criteria to be worse for late types. 



however our selection function also takes into ac- 
count the availability of rotation width data which 
selects against earlier types spirals. Without a re- 
liable luminosity function (LF) separated by mor- 
phological types along the spiral sequence there is 
no way to tell if the change in slope is due to incom- 
pleteness, or telling us something physical about 
the different types of spirals. There is also a possi- 
bility that it is related to variations in the amount 
of internal extinction in different types of spirals 
as a function of their intrinsi c size, althoug h that 
effect is not large in the NIR (jMasters et a/.l l2003'). 
In any case the difference in the relations between 
spiral types should be corrected for to construct 
the final "universal" TF relation for Scs. 

6. Cluster peculiar velocity. To combine all the indi- 
vidual cluster samples into a global TF relation, a 
peculiar velocity must be estimated for each clus- 
ter. We use the peculiar velocities derived in M06 
which considered the same set of clusters as is used 
here. For more details on this step we refer the 
reader to M06. Below we discuss further the im- 
pact this has on the zeropoint of the TF relation, 
and possible error introduced into the zeropoint by 
using a subset of the most distant clusters to define 
a zero-point rest frame. 

The global template TF relation is shown in Figure 
[2] for the in+ samples in all three bands. A summary 
of linear fits (direct, inverse and bivariate respectively 
considering errors in the magnitude direction, width di- 
rection and in both directions) to various subsets of the 
in+ sample is given in Table [21 In that table zeropoints 
(labeled a) are reported for Ho — 100 km s~^ Mpc~^. 
The slopes are labeled, h. The last four columns of Ta- 
ble [21 list the error on the bivariate zeropoint and slope, 
and the standard deviation and absolute deviation of the 
points from the bivariate relation. We list results for 
three fitting methods to aid comparisons with other re- 
sults {e.g. Section [3|), but our preferred method is the bi- 
variate fit which considers errors in both the magnitude 
and rotational width. Direct fits consider errors only in 
the magnitude and are the most affected by magnitude 
incompleteness in the sample. It is often suggested that 
using inverse fits completely removes the effects of in- 
completeness bias - however since inclination dependent 
corrections for both the magnitudes and widths create 
implicit cross-dependences in the final quantities there 
will be some incompleteness bias still present in the in- 
verse fits. Our preferred result (using the bivariate fit 
and all corrections described above) for the TF relation 
is then: 

-51og/i = -22.030 - 10.017(logVF-2.5), 
Mh - 5 log = -21.833 - 9.016(logM^ - 2.5), 
Mj - 5 log /i -20.999 - 9.070(log W - 2.5). 
In the fits both the observational errors and an esti- 
mate of the intrinsic error of the TF relation are used to 
weight points. We do not fit for the intrinsic scatter as 
the same time as the slope and zero-point, but measure 
it after and iterate the fit and scatter measurement until 
convergence is reached. 

The zeropoint derived for the TF relation from the 
"basket of clusters" method comes from an assumption 
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TABLE 1 

K-BAND Incompleteness Function Co-efficients and Cluster Size Bias. 
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Redshift 
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mags 






mags 
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N383 


4865 


14 


42 


-21.03 


1.01 


0.03 


-0.002 


-0.010 


N507 


4808 


17 


28 


-21.87 


0.42 


0.03 


-0.002 


-0.010 


A262 


4665 


28 


65 


-22.56 


0.69 


0.03 


-0.003 


-0.010 


A397 


9594 


13 


13 


-23.19 


0.49 


0.02 


-0.001 


-0.001 
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A569 


6011 


13 


15 


-20.43 


0.38 


0.02 


-0.001 


-0.006 


A634 


7922 


6 


6 


-21.65 


0.49 


0.02 


-0.001 


-0.001 


Cancer 


4939 


22 


75 


-21.88 


0.66 


0.03 


-0.002 


-0.009 


A 779 


7211 


13 


19 


-22.85 


0.47 


0.02 


-0.001 


-0.004 


A1314 


9970 


6 


6 


-20.70 


0.05 


0.02 


-0.001 


-0.001 


A1367 


6735 


37 


54 


-23.41 


0.62 


0.02 


-0.001 


-0.005 


Ursa Major 


1101 


16 


25 


-17.26 


1.24 


0.14 


-0.046 


-0.186 


Coma 


7185 


36 


52 


-22.59 


0.42 


0.02 


-0.001 


-0.004 


A2iy9/7 


8996 


8 


25 


-23.18 


0.43 


0.02 


-0.001 


-0.003 


Pegasus 


3519 


20 


39 


-19.73 


0.12 


0.04 


-0.005 


-0.018 


A2634 


8895 


18 


35 


-22.48 


0.41 


0.02 


-0.001 


-0.003 




7867 


9 


9 


-23.07 


0.31 


0.02 


-0.001 


-0.001 


Ai877 


6974 


7 


7 


-21.51 


0.36 


0.02 


-0.001 


-0.001 


A 194 


5036 


24 


32 


-21.24 


1.14 


0.03 


-0.002 


-0.009 


Eridanus 


loot) 


io 


z / 


-io. ( O 


U.oo 


U.IU 


n no/1 


n nnc; 


Fornax 


1321 


11 


23 


-17.58 


0.43 


0.11 


-0.032 


-0.129 


A496 


9809 


4 


4 


-22.66 


0.80 


0.02 


-0.001 


-0.001 


Antlia 


3120 


12 


29 


-21.61 


0.18 


0.05 


-0.006 


-0.023 


Hydra 


4075 


23 


32 


-21.61 


0.45 


0.04 


-0.003 


-0.014 


NGC 3557 


3120 


5 


11 


-21.79 


0.62 


0.05 


-0.005 


-0.020 


Cen 30 


3322 


20 


33 


-21.37 


0.80 


0.05 


-0.005 


-0.020 


ESQ 508 


3196 


7 


15 


-21.88 


0.65 


0.05 


-0.006 


-0.022 


A3574 


4817 


10 


25 


-22.40 


0.77 


0.03 


-0.002 


-0.010 


Pavo II 


4444 


7 


24 


-23.92 


0.31 


0.03 


-0.003 


-0.011 


Pavo I 


4055 


6 


14 


-21.06 


0.21 


0.04 


-0.003 


-0.014 


MDL59 


2317 


7 


15 


-20.06 


0.39 


0.07 


-0.010 


-0.042 



Note. — Column(l) - cluster name ordered by RA in the North (top half) and South 
(bottom half). Column(2) - cluster redshift in the CMB frame. Columns (3) & (4) - 
number of galaxies in in and in-|- samples respectively. Columns (4) & (5) - parameters 
of the incompleteness function for the cluster (Equation |2j • Column (6) - estimate of 
size cTd, divided by an estimate of the distance, d = Vcmb/Ho, for each cluster. Columns 
(7) & (8) - cluster size— incompleteness bias for the in and in-|- samples respectively. 

that the net motion of a subset of the most distant clus- 
ters is zero in the CMB frame (the so-called "cluster ref- 
erence frame" ) . We use in this paper the same zeropoint 
calibration (in terms of the peculiar velocities of the clus- 
ters in the CMB frame) that was derived in M06 for the 
I-band TF relation. This zeropoint could have a system- 
atic offset from the true zeropoint if there is a net motion 
of the cluster reference frame with respect to the CMB. 
The expected size of this possible offset was estimated 
in M06 under the assumption that the peculiar veloci- 
ties of clusters are characterized by an rms ID velocity 
dispersion with = 298 ± 34km (as measured in 
M06). They use this dispersion to estimate the possible 
offset from the random peculiar velocities of clusters to 
be I AM| 0.03 mag. This would be an underestimate of 
the zeropoint offset if there were coherent peculiar mo- 
tions on the scale of the entire cluster sample. How- 
ever as the sample extends to a distance of cz = 10000 
km s""*^ and clusters are located all over the sky a large 
effect from coherent peculiar velocities seems unlikely. 
In M06 a Cepheid calibration of the I-band TF relation 
was compared to the "basket of clusters" zeropoint and 
agreement is found for a Hubble constant of iJo = 74 ± 2 
(random) ±6 (systematic) km s~^ Mpc~^. The fact that 
this is identical to H p = 74 ± 2 km s~^ Mpc^^ found by 
ISanchez et al. I (2006) from a combination of WMAP and 
the 2dF galaxy redshift survey suggests that the "basket 
of clusters" zeropoint is not biased significantly by any 



-24 
-16 



^-24 

isfl : 
0-22 - 

'«-20 - 

^-18 ^ 

"''-IB El 



(a) K band 




(b) H band 

■ • v?.:'.\ 




^-24 

-18 T 
S-16 ^- 



( c ) J band 



.• r . •• 




2.2 



2.4 
log W 



2.6 



2.8 



Fig. 2.— Global TF relation for (a) K-band, (b) H-band (c) J- 
band(in-|- sample). The solid line is all panels shows the respective 
bivariate fit to the data (see Table 2). 

bulk flows across the sample. 

3.1. Wavelength Dependence 

As found by numerous other authors (see Section 5 
below for a discussion of some recent results in the NIR) 
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TABLE 2 

TF PARAMETERS 



Sample N 



Direct fit 



Inverse fit 

^inv bin 



Bivaxiate fit 



O'SD O-abs 



K-band 



No corrections 



Full 888 


-22.271 


-6.548 


-22.261 


-10.032 


-22.215 


-7.254 


0.019 


0. 


110 


0.571 


0.401 


Morphological correction only 




















Full 888 


-22.312 


-7.932 


-22.250 


-10.005 


-22.228 


-8.919 


0.017 


0. 


102 


0.510 


0.349 


Bias correction only 






















Sa 101 


-22.461 


-5.473 


-22.297 


-7.565 


-22.143 


-6.954 


0.048 


0. 


239 


0.463 


0.337 


Sb 345 


-22.132 


-7.205 


-22.227 


-9.417 


-22.034 


-8.635 


0.025 


0. 


165 


0.513 


0.353 


Sc 442 


-22.031 


-8.329 


-22.409 


-11.530 


-22.039 


-10.092 


0.022 


0. 


,175 


0.579 


0.427 


Full 888 


-22.094 


-7.296 


-22.252 


-10.011 


-21.997 


-8.647 


0.015 


0. 


093 


0.551 


0.398 


All corrections 
























Full 888 


-22.107 


-8.849 


-22.244 


-9.996 


-22.030 


-10.017 


0.015 


0. 


098 


0.523 


0.372 


Low mass 303 


-21.831 


-8.280 




















High mass 374 


-22.281 


-8.058 




















H-band 


No corrections 
























Full 888 


-21.990 


-6.438 


-21.981 


-9.595 


-21.943 


-7.068 


0.019 


0. 


,117 


0.580 


0.414 


Morphological correction only 




















Full 888 


-22.012 


-7.608 


-21.972 


-9.580 


-21.954 


-8.289 


0.018 


0. 


114 


0.542 


0.382 


Bias correction only 






















Sa 101 


-22.196 


-5.383 


-22.004 


-7.429 


-21.996 


-6.084 


0.054 


0. 


265 


0.475 


0.355 


Sb 345 


-21.910 


-6.833 


-21.954 


-9.100 


-21.842 


-7.803 


0.027 


0. 


165 


0.526 


0.365 


Sc 442 


-21.838 


-7.807 


-22.119 


-10.976 


-21.847 


-9.165 


0.024 


0. 


,189 


0.569 


0.423 


Full 888 


-21.879 


-7.001 


-21.975 


-9.587 


-21.812 


-7.896 


0.016 


0. 


105 


0.553 


0.401 


All corrections 
























Full 888 


-21.893 


-8.235 


-21.967 


-9.575 


-21.833 


-9.016 


0.016 


0. 


,106 


0.530 


0.382 


Low mass 303 


-21.612 


-6.993 




















High mass 374 


-21.972 


-7.838 




















J-band 


No corrections 
























Full 888 


-21.236 


-6.343 


-21.218 


-9.448 


-21.192 


-6.905 


0.020 


0. 


133 


0.613 


0.446 


Morphological correction only 




















Full 888 


-21.252 


-7.433 


-21.208 


-9.427 


-21.198 


-8.030 


0.019 


0. 


,120 


0.587 


0.423 


Bias correction only 






















Sa 101 


-21.349 


-5.604 


-21.234 


-7.347 


-21.218 


-6.093 


0.056 


0. 


297 


0.494 


0.363 


Sb 345 


-21.064 


-7.084 


-21.188 


-8.914 


-21.016 


-7.795 


0.027 


0. 


,185 


0.562 


0.404 


Sc 442 


-21.021 


-8.244 


-21.352 


-10.831 


-21.017 


-9.228 


0.026 


0. 


,202 


0.625 


0.476 


Full 888 


-21.037 


-7.408 


-21.207 


-9.425 


-20.985 


-8.104 


0.017 


0. 


,117 


0.597 


0.444 


All corrections 
























Full 888 


-21.050 


-8.428 


-21.199 


-9.411 


-20.999 


-9.070 


0.017 


0. 


,116 


0.578 


0.426 


Low mass 303 


-20.695 


-6.588 




















High mass 374 


-21.086 


-8.289 





















Note. — Column (1) describes the sample used to fit to, column (2) shows the total number of galaxies in that 
sample. Columns (3-8) show parameters for M — 5 log /i — a + b{logW — 2.5) using direct, inverse and bivariatc 
fits. Also shown for the bivariatc fits arc the error on a and b (columns 9 & 10), and the scatter calculated as the 
standard deviation (column 11) and absolute deviation (column 12). Inverse and bivariatc fits arc not reported for 
the sample divided by galaxy size since the artificial log W cut-off which is applied creates a large bias in those 
fits. 
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Fig. 3. — Trends of the (a) zeropoint and (b) slope of the TF 
relation with wavelength. I-band results from the related sample 
of M06. Circles shows the result for Sas, squares for Sbs, triangles 
for Scs and stars for the combined sample (with magnitudes of Sas 
and Sbs corrected to fall on the Sc relation). 
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Fig. 4. — K-band TF relation for the in-|- sample separated by 
morphological type (a) Sa and earlier; (b) Sb; (c) Sc and later. 
In each panel open symbols show the magnitudes corrected for all 
biases except morphological type. The dashed lines show the fits 
to these points. The filled symbols are the magnitudes after the 
morphological type correction is applied and the solid line in all 
panels shows the final fit to the full sample. 

the TF relation steepens as the wavelength of the band it 
is measured in increases. We find no significant evidence 
for wavelength dependence of the scatter (see Section 4 
below). Figure [3] summarizes the trends of the slope and 
zeropoint with wavelength including both the relations 
reported in this paper for J, H and K-bands as well as 
the 1-band result from the related sample of M06. 

We find that the slope of the TF relation is very close to 
the simple theoretical expectation of & = —10 (L oc v^) 
in K-band, with Lk oc 7;4 00±o.o4 H-band we find 
Lu oc ^;3-6i±o.o4 j-band we find Lj oc u3-63±o.05_ 

For their related sample, and using the same method as 
in this paper M06 find that in I-band L/ oc ^;3-i4±o.o4_ 

The zeropoint of the TF relation also changes with 
band. This change as expected is consistent with the 
average colors found for galaxies in the template sample 
oiJ-K= 1.04 andH-K ^ 0.28, with uj-qk = 1.03 
and an — olk = 0.20. 

3.2. Morphological Dependence 

We find in all three 2MASS bands that the slope of the 
TF relation steepens as the sample changes from earlier 
type to later type spirals (see Figures [3] and [4]) . This is 



also as found in I-band for the related SFI-I--I- template 
sample (M06). The amount the slope changes appears 
to be relatively independent of the band of the TF rela- 
tion magnitudes. We find little evidence for a zeropoint 
change between Sb and Sc galaxies, but Sa's do appear 
to have a zeropoint 0.2 mags brighter than later type 
spirals in all three 2MASS bands. 

These findings are contrary to the commonly held idea 
that morphological variation of the TF relation disap- 
pear almost completely in the NIR {e.g. Noordermeer & 
Verheijin 2007). This conclusion w as based on a lack of 
morph ologi cal dependence found bylAaronson fc Mould I 
(flOO and lPeletier fc Willned (flQOl in the H-band TF 
relation, and was used to argue that the morphological 
dependence seen in the bluer bands could be attributed 
to differences in star formation history between types of 
spirals. The clear (if small) morphological dependence 
we pick up is presumably a result of our much larger 
sample size, and indicates that the change in the TF re- 
lation cannot be totally accounted for by the details of 
star formation history and instead may be related either 
to something more fundamental about the dynamics of 
disk galaxies along the Hubble sequence, or possibly to 
the dust content or metallicity of different types of spi- 
ral galaxies (although dust content is not seen to have a 
large impact in the NIR, Masters et al. 2003) . In M06 we 
interpreted the change in slope between different types 
of spirals (if real and not an effect of varying incomplete- 
ness for different spiral types) as possibly being due to the 
dark matter halos of earlier type spirals being less con- 
centrated at a given rotation velocity than later types. 
The large caveat here is that if the incompleteness bias 
is significantly larger for Sas than later type spirals, the 
slope of the Sa relation may be artificially shallowed and 
the zeropoint artificially brightened by the effects of this 
bias. Naively, in the NIR we expect the selection criteria 
to be worse for late types, however our selection function 
also takes into account the availability of rotation width 
data which selects against earlier types spirals, so either 
could be true. However it should also be noted that much 
of the change in slope is driven by massive Sa's (and Sb's) 
which should be unaffected by the incompleteness bias. 

3.3. Luminosity Dependence 

There is present in Figure [2] a hint that both at the 
very high mass and very low mass ends of the TF relation 
galaxies turn off slightly from the mean trend. We divide 
the full sample into "luminosity classes" to study this 
further using logW^ limits of logW^ < 2.4; logW^ > 2.5 
which roughly selects the lowest and highest mass thirds 
of the sample. We fit the TF relation to these subsets; 
the results are in Table [2l In this case we are adding in 
an artificial strict rotation width limit so biasing both 
the inverse and bivariate fits. We therefore report only 
the direct fit in TableO We see evidence that at both the 
high and low mass ends of the relation there is a slightly 
shallower slope, however in both cases and in all 3 bands 
this difference is very small. 

Disentangling fully the luminosity dependence of the 
slope and the morphological dependence will be difficult 
since as illustrated in Figure |4] the log W distribution of 
the different types of spirals is quite different with many 
more high mass Sas compared to Scs. We note here that 
we look for luminosity dependence after correcting for 
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morphological type as described in Section 3. Searching 
for a variation is also complicated by the impact of the 
incompleteness bias correction. We assume a linear TF 
relation to calculate the bias, which has a large effect at 
the low mass end of the relation. Separating a variation 
of the TF at the low mass from this incompleteness bias 
is obviously not trivial. 

There have been several studies in the optical pass- 
bands of the low mass end of the TF relation including 
Kannappan, Fabricant & Franx (2002) who describe a 
split in the dwarf extension of the TF relation in U, B 
and R-bands confirming findings by Pierini (1999) in the 
B-band. Hints of a deviation can also be seen in the 
Tband TF presented in M06. Typically the deviation 
has been seen to make the slope of the linear relation 
steeper for low mass galaxies, while here we see hints of 
it shallowing. However in other studies (which may have 
been more complete for low mass galaxies) there was no 
correction for incompleteness bias, which is found to be 
quite large in our sample. The bias correction we apply 
brightens all points at log W < 2.4 by up to a magnitude 
or more (see Figure [1]), so any conclusions we draw on 
the low mass end of the relation are highly dependent 
on the assumption which went into the bias correction. 
The choice of rotation width measure and turbulent mo- 
tion corrections applied to the widths also impact the 
small rotation width galaxies much more greatly. We 
use rotation widths based on the width of the HI line at 
50% of the peak, and corrected for a turbulent motion 
using a linear correction of 6.5 km as justified in 
ISpringob etoT (2005). 

At the high mass e nd o f the TF relation, 
iNoordermeer fc Verheiien] (|2007D recently argue that 
very massive galaxies lie systematically to the right of 
the mean K-band TF relation. This the same as the 
trend we see here and which is also present in the I-band 
relation of M06. 

4. THE SCATTER IN THE TF RELATION 

The scatter in the TF relation is important, both to 
understand the errors on distances derived from TF, and 
as a constraint for models of galaxy formation. As in 
G97b and M06 we observe that the scatter in the TF 
relation increases with decreasing rotation width. We 
therefore parametrize the scatter as 

a ^ a + b(\ogW -2.5). (3) 

In both G97b and M06 there is a discussion of the effect 
incompleteness bias can have on the observed scatter. We 
find that in this sample the bias on the scatter is fairly 
small, so we neglect to correct for it. It should then be 
noted that the scatters quoted here are underestimates 
of the real scatter by something on the order of 10% 
at the low mass end of the relation. At the high mass 
end (above logW^ ~ 2.5 or W ^ 320 km s~^) there is 
no longer any impact from incompleteness bias on the 
scatter. 

The measured scatter (standard deviation) is shown 
for all 3 2MASS bands in Figure [5l A linear fit to the 
total scatter gives 

eobs,K = 0.54 - 0.87(log 1^-2.5) 
eobs3 = 0.54 -0.89(logW^- 2.5) 
eobs,j = 0.56 - 0.66(logW^- 2.5) 



(shown by the dash-dot lines in the figure) . We estimate 
the total measurement errors from the template sample 
as a function of rotation width. These values are shown 
by the dotted lines in Figure[5l The lower line at a ~ 0.1 
shows the measurement error on the total magnitudes (in 
the respective bands) which at all widths is the least im- 
portant source of error. The middle line shows a value of 
|6|(JvK which is the error on the measured rotation widths 
expressed in magnitudes. Adding these two contribu- 
tions in quadrature approximates the total measurement 
error (this neglects the covariances between the two val- 
ues which arise from the inclination corrections). Using 
this estimate of the total measurement error we also fit 
for an intrinsic scatter which varies with rotational width 
finding: 

ei„t,K= 0.35 -1.17(logiy- 2.5) 
ei„t H = 0.38 - 1.14(log W - 2.5) 
ei„t J = 0.39 - 0.78(log - 2.5) 

These estimates of the intrinsic scatter in the TF rela- 
tion are shown by the dashed line in Figure [Sj the solid 
line shows the quadrature sum of this and the total mea- 
surement error. This sum provides a better fit to the 
observed errors than the straight line fit. 

There is no significant evidence for a variation of the 
TF scatter with band, however there is a consistent in- 
crease in the zeropoint and decrease in the variation with 
logVF from K- to J-band. We also find a larger total 
scatter than in the I-band sample of M06 (which reports 
ei„t,i = 0.35 - 0.37(logVF - 2.5)) but this may be re- 
lated to the larger errors on correcting the widths for 
inclination by using 2MASS axial ratios. This may not 
be totally accounted for in our estimate of the total mea- 
surement error and certainly could add significant covari- 
ance between the magnitude and width errors, so a direct 
comparison between the scatter measured here and that 
of M06 should be done carefully. 

In the sub-samples divided by morphological type we 
have not fit separately for the intrinsic scatter a a func- 
tion of log W. From Tabled however we can compare the 
total scatter of the different subsets from their own TF 
relations. We find in all three bands a consistent increase 
of the total scatter from Sa to Sc galaxies. However dis- 
entangling this effect from the luminosity dependence of 
the TF scatter will be tricky. There are many more low 
mass Scs than Sas and we measure more intrinsic scatter 
at the low mass end of the TF relation, so the average Sc 
should have a larger deviation from the mean TF relation 
than the average Sa. Whether the causal reason for this 
larger scatter at the low mass end of the TF is the lower 
mass of the galaxies, or the increased number of later 
types is beyond the scope of this paper, but provides an 
interesting question for further study. 

5. COMPARISON WITH OTHER TF RELATIONS 

The first calibration of the TF re lation in NIR bands 
was done bv lAaronson et aC] ()1979( ). They found a slope 
of approximately b = —9.5 (or L oc v^'^) in the H- 
band using samples of galaxies in the two nearby clus- 
ters, Virgo and Ursa Major. This result was similar to 
the naive expectation of L oc u"* which arises from a 
simple consideration of the dynamics of spiral disks and 
the assumption that M/ L is constant (see Appendix C 
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Fig. 5. — Total observed scatter in the TF relations in (a) K-band, (b) H-band and (c) J-band. The dash-dot line shows a linear fit to 
the total scatter. Also plotted is the error budget for all galaxies in the dataset as a function of rotation width. The dotted lines show 
errors associated with the photometry and width measurement (multiplied by the TF slope to be expressed in magnitudes). This is added 
in quadrature to a width dependent intrinsic scatter (shown by the dashed line — the sum is shown by the solid line) which is fit for. 

lations uncorrected for incompleteness bias should all be 
shallower than what we derive here (with no other differ- 
ences in the sample make up) - how much shallower they 
are will depend on the degree of incompleteness in the 
sample. Then there are hints that the TF relation slope 
varies with galaxy luminosity meaning that samples with 
different distribution of galaxy luminosities may also give 
different TF relations. With all these different factors 
to worry about it's clearly over simplistic to compare 
different calibrations directly. Never-the-less we provide 
such a direct comparison in Table[3l with brief comments 
about the different sample make-ups. In the remainder 
o f this section we d iscuss each calibration in more detail. 

iTullv fc PieTcel (j200d!) derive TF relations in the B, 
R, I-bands using samples of galaxies drawn from 5 
nearby clusters (with distance c alibrators for th e zero - 
point). In a companion paper iRothberg etHl (|2000r ) 
used a similar sample to derive the K-band TF rela- 
tion. These papers claim that using the inverse re- 



of ICourteau et al.ll2007l for a recent dis cussion of simple 
analyt ic predictions for the TF slope). lAaronson et al. I 
(|1979|) also observed that the scatter was lower than in 
the B-band TF, paving the way for significant work using 
the NIR TF relation. 

There has been quite a lot of work on deriving template 
TF relations in the NIR bands since this initial paper. 
Here we attempt to compare some of these derivations 
with what we find here. This comparison is not complete, 
and is biased towards recent calibrations of the NIR TF 
relations. As has been shown here (and in M06) both the 
slope and zeropoint of the TF relation vary with morpho- 
logical types. If this dependence is not corrected for then 
samples with different morphological type distributions 
should show different TF relations. In particular cluster 
samples may be more dominated by earlier type spirals 
(so expected to have shallower TF slopes) than field sam- 
ples. Also it is not yet common practice to correct for 
incompleteness bias when deriving TF relations. TF re- 
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TABLE 3 

COMPARISION OF SOME PUBLISHED NIR TF RELATIONS. 



Reference 


I-band 


J-band 


H-band 


K-band 




Comments 




a 


b 


a 


b 


a 


b 


a 


b 






This work 






-21.65 


9.07 


-22.49 


-9.02 


-22.68 


-10.02 


N = 


890 (bias corrected, Sc) 


This work 






-21.67 


7.80 


-22.50 


-7.80 


-22.69 


-8.63 


N = 


345 (bias corrected, Sb) 


This work 






-21.87 


6.09 


-22.65 


-6.08 


-22.80 


-6.95 


N = 


101 (bias corrected, Sa) 


(1) 


-21 47 


-8.59 










-21.99 


-9.29 


Ni = 


= 1303, Nk = 360 








-22.1 


-6.3 


-22.6 


-6.4 


-22.9 


-6.6 


AT 


iuuu unDiabeQ , o^/o oaj 


i 


-21.50 


-7.85 














N = 


807 (bias corrected, Sc) 




-21.38 


-7.07 














N = 


281 (bias corrected, Sb) 




-21.42 


-5.52 














N = 


61 (bias corrected, Sa) 






-6.37 




8.20 




-8.71 




-9.02 


N = 


436 (edge-on, field) 


g! 


-21.51 


-9.6 










-23.05 


-10.6 


N = 


16 (UMaj, complete?) 












-22.01 


-10.5 


-22.34 


-10.4 


N = 


19 (Cepheids) 




-21.57 


-8.11 














N = 


155 (cluster) 
















-23.17 


-8.78 


N = 


65 (82 % Sb/Sc) 




-21.01 


-7.68 














N = 


555 (bias corrected, Sc) 


(10) 












-9.5 






N = 


29 


Note. — 


Previous 


calibrations of the 


TF 


■elation in 


I, J, H 


and K-bands are 


shown 


in order of publication date 



Zeropoints are adjusted for — 74 km 
tion from primary distan ces. References ~ 

IKarachentsev et al.l | [2q0^. (5) -IVerheiienI II2001I) . (6) - iMacri 12001,). (7) 
(120001 ). (9) G97b. (10) Uaronson ef aL I I I1979I) 



M pc where necessary for comparison ^ 
- ICour teau et al. (2007), (2) Theuroau et 
- TuUy & Pierce 



vitli those based on calibra- 
al. (2003), ( 3) - M06, (4) - 
(2000), (8) - IRothberg et al.l 



lation "nulls" any bias introduced by sample incom- 
pleteness (however as discussed in several places this is 
only true in the case that there is no covariance be- 
tween luminosity and width errors which is difficult to 
achieve in practice because of inclination dependent cor- 
rections) They find in I-band (for 155 galaxies) Mj = 
-21.57 - 8.11(logM^ - 2.5), while in K-band [N = 65) 
they find = -23.17 - 8.78(log IF - 2.5). This I-band 
relation is steeper than the bias corrected relation of M06 
{N = 807), although consistent with M06's inverse fit 
with no incompleteness bias cor rection of Mj — 5 log h = 
-20.89 - 8.341(logir - 2.5). IRothberg et "all (l2000)'s 
K-band relation is much shallower than both our pre- 
ferred K-band result and our inverse fit to the sam- 
ple with no incompleteness bias correction applied of 
Mk - 51og/i = -22.250 - 10.005(logVF - 2.5), which 
could be due to the relatively small K-band sample size 
or the effects of incompleteness. This relation is similar 
to our fit to Sb galaxies, so the differe nce could also be 
explained if the IRothberg et al.l (|2000( ) sample is domi- 
nated by Sb and earlier type spirals. However this is not 
the case - this sample actually has almost equal numbers 
o f Sb and S c galax ies (see their Table 1). 

IVerheiienI (|2001h derive B, R, I and K-band templates 
for a sample of galaxies in the Ursa Major cluster. They 
claim this sample {N = 16) is complete. They find 
Mi = -21.51 ± 0.83 - (9.6 ± OA){logW - 2.5) and 
Mk = -23.05 ±0.95- (10.6 ±0.4)(logVl^- 2.5). The K- 
band relation is very similar to that found in this paper. 
The I-band relation is much steeper than that found in 
M06 and G97b. 

Macri (2001) provides a Cepheid calibration of the 
H and K-band TF templates using 19 galaxies with 
Cepheid distances, for the 20th magnitude per arcsec- 
ond squared isophotal magnitudes in 2MASS they find 
Mk = -22.34 ± 0.07 - 10.4 ± 0.5(logW^ - 2.5) and 
Mh = -22.01 ± 0.02- 10.5 ±0.4(logiy- 2.5). As must 
be the case for a sample based on galaxies with primary 
distances the sample size is small, and since galaxies with 
Cepheid distances tend to lie preferentially at the high 
mass end of the TF relation usually it's not a good idea 



to derive a slope from such a sample. The Macri (2001) 
K-band slope agrees with our result quite well, while the 
H -band slope is somewhat s teeper. 

IKarachentsev et al.l ()2002| ) derive B, I, J, H and K TF 
relations for the most edge-on spirals in 2MASS. The 436 
galaxies used are field spirals (therefore likely dominated 
by late type spirals) and are placed on the TF relation 
using absolute magnitudes derived with no correction for 
peculiar velocities. This could bias the slope as the im- 
pact of peculiar velocities is likely to be larger for galax- 
ies at the low width end of the relation since they are 
more likely to be nearby. There is also no attempt to 
correct for sample selection biases. Only the slopes are 
reported: in the NIR they find s/ = —6.37 ± 0.18, sj — 
-8.20±0.24, sjj = -8.71±0.25, sk = -9.02±0.25, shal- 
lower than we find in our bias corrected work as should 
b e expected. 

iTheureau et al.l ([2007^ use the "unbiased plateau" of 
field galaxies (as described in Theureau et al. 1998) to 
derive TF relations in the 2MASS J, H and K-bands. 
This method is iterative in that an input TF is used 
to derive distances to the galaxies from which an unbi- 
ased subset is used to construct a new input TF... and 
so on. They find (converted to the terminology used 
in this paper), Mk = -23.5 - 6.6(logT4^ - 2.5), with 
cr = 0.45; TV = 990, Mh ^ -23.2 -6.4(log 1^-2.5), with 
a = OAT, TV = 1166 and Mj = -22.7 - 6.3(log W - 2.5), 
with a = 0.46, N = 960. A Hubble's constant of Ho = 57 
km Mpc~^ 
rect to Ho = 74 km s~ 
other work (see Table [3]). The slopes of these relations 
are significantly shallower than we find in this paper 
for our Sc sample (and that are reported elsewhere in 
the literature). A possible explanation is that the dif- 
ference is related to the lumino sity distribution of the 
sample, in that the sample of ITheureau et al.l (|2007l ) 
acts like a volume limited sample and is therefore dom- 
inated b y massiv e spirals , while the s amples in this pa- 
per. in iVerheiien. (i200l1). ITuIIv fc Pierc e (2000) and in 
IKarachentsev et"al1 ( 20021 ) have many more low mass 
galaxies in them. If lower mass galaxies exhibit a steeper 



is implicit in these zeropoints. We cor- 
^ Mpc~^ for comparison with 
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TF relation than massive spirals this could explain the 
difference between the measured slopes (Theureau 2007; 
priv. comm.). This idea is backed up to some extent 
by our investigation of the luminosity dependence of the 
slope, however these slopes are still significantly shal- 
lower than the direct fits to the most massive third of our 
sample suggesting that this cannot be t he full explana- 
tion fo r the difference. Also Figure 9 of iTheureau et ahl 
(|2007l ) shows that the sample to which the TF relation is 
fit in that paper does not differ significant ly in log W ex- 
tent from our sample (note that Theure au et al.ll2007l use 
log V = log W — log 2 ^ log W — 0.3), so it does not seem 
likely that this can explain the large difference in TF 
slope. What appears to be a more likely explanation for 
the difference in the slopes is the morphological m ake-up 
of the sample. A field sample like that in Theu reau et al.l 
(|2007l ) might be expected to be dominated by the latest 
type spirals, however the selection of only the brightest 
field gal axies should instead bi as it to earlier types. In 
fact the ' Theureau et al.l (|2007[ ) sample (or at least the 
galaxies with new HI observations listed in their Table 
3) is composed of more than half (64 %) Sab and ear- 
fier types (T < 2), and only 19% Sbs, and 17% Scs and 
later. Comparing the slopes they find with our "Sa" TF 
relation the agreement becomes reasonably good. 

Courtcau et al. (2007) derive I- and K-band TF rela- 
tions from a large dataset of both field and cluster spirals 
constructed from the combination of 4 smaller samples. 
In total they have 1303 galaxies which are placed on the T 
band TF relation, only 360 of which ("the brightest SCII 
galaxies") are used to derive the K-band relation. They 
find (converted to the terminology used here, and to 
Ho = 74kms-i Mpc-i) Mj = -21.47-8. 59(log 1^-2.5) 
and Mk = -21.99 - 9.29(logT4^ - 2.5), fairly similar to 
the relations in M06 (for I-band) and found here (for 
K band). They also notice trends in the slope and ze- 
ropoint with both band and morphological type of the 
sample of the same sense we find in this paper (and was 
found in M06 for I-band) but of a smaller size. The 
slope of their Sc Tband TF relation {N = 505) is -9.19 
mags/dex while for Sa's {N =117) they find a slope of 
-8.25 mag/dex. In K-band they find for Sc's {N = 166) 
a slope of -9.80 mags/dex and for Sa's {N = 56) a slope 
of -9.12 mags/dex. 

It is interesting to ask if our own Milky Way galaxy 
lies on th e TF relation we find i n this paper. A previous 
study by ' Malhotra et al.l (|i996') showed that the Milky 
Way was within 1.5cr of the NIR TF relations of other 
nearby galaxies. The availability of photometrically ac- 
curate all-sky surveys in these bands (like 2MASS and 
DIRBE) means that it is easier (although still difficult) 
to estimate the total luminosity of the Milky Way in 
these bands to place it on the TF relation. The total 
luminosity of the Milky Way w as estimated in t he J a nd 
K bands from DIRBE data bv lMalhotra et all ()1996f ) as 
Mk = -24.06 mags and Mj = -23.05 mags. Recent 
estimates of the total rotation veloci ty of the Milky Way 
vary between v = 226 ± 15km s^-^dLepine et al.|[2008f) 
and V = 270km s~^ ()Mendez et al.lll999D . We therefore 
use W = 496 ± 44 km s-\ or \ogW = 2.69 ± 0.05 to 
cover both ends of the estimate. The best guess of the 
morphologic al typ e of the Milky Way is that it is an Sc or 
SBc (.Hodg^ ligsl . meaning that it should lie on the Sc 
relations we derive above (see Equations 4) with no cor- 



rection for morphology. We find that for Hq — 74 km s^ 
Mpc^^, the Milky Way is slightly under luminous, but 
consistent within the measured intrinsic scatter of both 
the J and K-band TF relations, being 0.32 mags dimmer 
than the J-band TF relation, and 0.53 mags dimmer than 
the K-band TF relation. If the larger rotational speed is 
accurate the deviation from the TF relations would be 
greater, using the small end of the rotation speed esti- 
mates the Milky Way is entirely consistent with lying on 
the TF relation derived here. 

6. SUMMARY AND DISCUSSION 

The 2MASS TuUy-Fisher Survey (2MTF) aims to mea- 
sure TuUy-Fisher (TF) distances to all bright inclined 
spirals in the 2MASS Rcdshift Survey (2MRS). Essen- 
tial to this project is a universal calibration of the TF 
relation in the 2MASS J (1.2 /^m), H (1.6 ^m) and K 
(2.2 //m) bands. We present the first bias corrected (or 
universal) TF template in these bands. Previously pub- 
lished NIR TF relations have either ignored the issue of 
incompleteness bias affecting the slope and zeropoint de- 
rived, or have attempted to construct a complete sample 
from a cluster, or by using only part of their data. 

As in earlier work, we find that the slope of the TF re- 
lation becomes steeper as the wavelength increases. We 
find that the slope is close to the simple theoretical expec- 
tation of L cx in K-bands, while is is slightly shallower 
at H- and J-band at L oc v'^'^. The change in zeropoint 
between bands in consistent with the mean colors of the 
galaxies in the template sample. 

We investigate the dependence of the TF relation on 
galaxy morphology showing that in all three bands the 
relation is steeper for later type spirals which also have 
a dimmer TF zeropoint. We correct for this dependence 
in deriving a global TF relation for Sc galaxies in the 
NIR. We also investigate the luminosity dependence of 
the slope of the TF relation, finding hints that at both 
the highest and lowest masses the slope of the TF relation 
is shallower than the mean relation - however this effect 
is not large. 

We fit the TF relation using direct, inverse and bi- 
variate fits. We favor the bivariate fit which takes into 
account errors on both the width and magnitude mea- 
surements (we include the measured intrinsic scatter in 
the weighting of the fits) finding a TF template for Sc 
galaxies (but which can be applied to all spirals with a 
morphological correction) of: 



Mk 


— 5 log h = 


-22.030- 


10.017(logM^- 


-2.5), 


Mh 


— 5 log h = 


-21.833 - 


9.016(logiy- 


2.5), 


Mj 


— 5 log h = 


-20.999 - 


9.070(logVF- 


2.5). 



These relations are bases on 2MASS extrapolated mag- 
nitudes and rotation widths measured from the HI line 
of galaxies at the 50% level of the (peak— rms) on a fit to 
the edges of the profile (Wfso; supplemental widths from 
optical rotation curves are corrected to the W^fso scale 
using the prescription described in M06 based on work 
of CatincUa et al. 2007). To measure a TF distance an 
observation of the inclincation, i, of the galaxy is also re- 
quired to correct the observed Doppler broadened width 
to the full rotational width of the galaxy, and to correct 
the observed total magnitude for extinction from dust 
internal to the galaxy. One would measure the distance 
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to a galaxy using the deviation of a galaxy from this 
TF relation (AM) and its measured recessional velocity 
using 

d = — 10 5 , (5) 



Ho 

or (equivalently, but now independent of Ho) the peculiar 
velocity in the CMB frame is 



Wpcc,CMB = WCMb(1 - 10 ). 



(6) 



Here AM — A/obs - M{W) is the difference between 
Mx{W) — 51og/i from the TF relations above (Equa- 
tions 4, where X = J, H or K describes the band the 
magnitude is measured in) and the corrected absolute 
magnitude of the galaxy, calculated from the observed 
apparent magnitude using 

Mcorr-51og/l = niohs-Ax-Ix+kx-Tx-^logVcMB-^^. 

(7) 

The observed magnitude of a galaxy above, has a cor- 
rection for extinction due to dust in our Galaxy (Ax), 
and internal to itself {Ix), as well as a cosmological K- 
correctio n (kx)- We use the prescription described in 
detail in Mast ers et al\ (2003) to apply these corrections 
to the 2MASS extrapolated magnitudes. Also needed is 
a morphological type correction, Tx{W) to correct for 
change in slope and zeropoint we see between different 
types of spirals. This correction depends on both the 
type, T and the rotational width of the galaxy and is 
described in Section 3 above. The observed width, Wfso 
must also be corrected for instrumental effects (Ag), in- 
clination, cosmological broadening and the effect of tur- 
bulent motions (At) in the disk: 



1 + z 



1 

sini ' 



(8) 



The instrumental corre ction will depend on t he details 
of the telescope used see lSpringob et al] (|200^ for more 
information on this. Based on arguments in Springob et 
al. (2005) we use At = 6.5km s"!. 



The scatter in the TF relation is important to under- 
stand, both to properly derive distance/peculiar velocity 
errors derived from TF, and as a constraint to models 
of galaxy formation. We fit for an intrinsic TF scat- 
ter which varies with galaxy luminosity (following the 
methodology of G97b and M06). We find no evidence for 
significant variation of the scatter with the wavelength 
of the magnitudes used. We do observe that later type 
spirals on the mean have a larger intrinsic scatter than 
earlier types, however disentangling this effect from the 
luminosity dependence of the scatter is tricky (since there 
are more later types at the low mass end of the relation 
where the scatter is observed to be larger) and beyond 
the scope of this paper. 

The TF relations derived in this paper will be used in 
future work to measure distances to galaxies in the 2MTF 
survey. The derivation of an unbiased or "universal" TF 
relation in the NIR J, H and K-bands should also be of 
significant use to constrain models of galaxy formation. 
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